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ABSTRACT 


This report reworks an earlier theory of the a 


early part of the radio flash from a low air burst 


by extending the results to quite general gamma ray 


vs. time histories, and to times about ten times 


longer than the earlier theory. 


The previous theory 


is also improved in that error estimates are given 


for all approximations used. 


are presented. 


I. INTRODUCTION 

Some time ago a method, called the high-frequen- 
cy approximation, was developed for calculating, ana- 
lytically, tne early part of the radio flash from an 


air burst. : Subsequently the method was extended 


to a ground parst.$*? 


The first of these references, 
(1) it 


is classified, and (2) in it many approximations were 


however, suffers from two defects, namely: 


made in order to obtain simple results without any 
estimates as to the limits of validity of the approx- 
imations. It has accordingly been deemed advisable 
to rework the theory removing as far as possible the 
defects. 

Our analytical theory of the radio flash re- 


quires three classes of approximations: 


1. Analytical approximations to the 
temporal and spatial behavior of 
the source functions, J and o. 

2. An approximate treatment of Maxwell's 


equations, arrived at by dropping 


certain troublesome terms. 


Typical signal shapes 


3. The evaluation of certain integrals 
which express the solution of the 


simplified set of Maxwell's equations. 


In Section II we derive the equations of the 
high~frequency approximation, first in full general- 
ity, and then specialized to the low-altitude air 
burst. In Section III we discuss the source func- 
tions J and o in sufficient detail to indicate tnat 
the analytic form we assume is sufficiently general 


Thus item 1. 


to cover all practical cases. above 
presents 


late the 


no problem. In Section III we also calcu- 
radial E-field which is needed in the next 
section. 

Section IV contains our results, namely ana- 
lytic expressions for the radiated signal, togetner 
with typical curves. These results are, of course, 
This 


approximation may be viewed as the first term of a 


based on our high-frequency approximation. 


series expansion of the solution to Maxwell's 
equations. We have made the usual heuristic esti- 
mate of the limitations of the theory by calculating 


the second term of this expansion and noting under 


what conditions it is negligibly small. The details 


of this calculation are in Appendix D. Generally 
speaking, we can be confident in the first micro- 
second of the calculated pulse. In order to express 
the solutions to our equations in simple closed form 
approximate expressions must be found for certain 
integrals. These approximations, together with re- 
mainder terms, are worked out in Appendices A and B. 
As the remainder terms given are exact, there can be 
no question as to the range of validity of our ap- 
proximations; at the expense of more complicated for- 
mulas, correction terms of any desired precision can 
be added to our formulae. 

Finally in Section V we discuss our results. 
In this section it is pointed out why one may ignore 
all nonlinear effects, as we have done throughout 


this work. 
Il. THE HIGH-FREQUENCY APPROXIMATION 


The Compton current, J, and the electrical 
conductivity, o, are both produced by a short pulse 
of ganma radiation whicn expands radially outward 
from the burst point with the speed of light. It is 
tnerefore appropriate to describe events in terms of 
proper time, +, rather than ordinary time, t 


(2.1) r= t-r/e, 


wnere r is distance from the burst point. When ex- 
pressed in terms of proper time, Tt, and position, r, 

~ 
the Compton current and the conductivity have approxi- 


mately the forms 


J= ir)f(r) 
(2.2) 1 
o = s(r)a(r) 


where f and g represent extremely short pulses. 


Crearly, 
A JJ ole 
(2.3) or™ -J/x and or o/x .. 


wnere A is the gamma-ray mean-free-path. Now a/c is 
about 2/3 microsecond and the pulses, f and g, are con- 
siderably narrower than this. Therefore the operator, 


J/or, operating on J or o, is much smaller than the 


operator, d/cdt. It is reasonable to assume that the 
same is true of the field quantities, at least for 
the early part of the signal. The high-frequency 
approximation consists in ignoring 3/9r of various 
field quantities compared with /oj1 of the same 
quantity. 


however, be done with circumspection, first noting 


The dropping of such small terns must, 


all cancellation that may exist among large terms 
before small ones are dropped. 

The easiest way to see how the high-frequency 
approximation works in detail is to write out the 
full set of Maxwell’s equations in polar coordinates, 
(r,0,), and in terms of the proper time, +, rather 
than ordinary time, t. Under this transformation, 


Maxwell’s equations take the form 


1 oF, 
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The first step is to eliminate By and Bo by formally 
integrating the last two equations 


=00 
= 
1 ) 
(2.5) ~ ‘sin 6 30 i E.cdr , 
“co 
T T 
rB. = rE, + oe E_edt - = E ed 
t ” 8 30 root ~ oF ere 
-@ -00 


Then substituting these results into the two 
equations for Ey and Eo gives 


0E 
r 


(2.6) 2 (rE,) + 2no( rE, ) = -enrd, + $ rs 


T T 
OB 2 2 
1 r 10 pO: Se a 
+B sin’ op ~ 2 droo Eedr + 5 dre rE car 
-@ Yee) 
and 
(2.7) 3 (rE_) + 2no(rE.) = -2nrJ_ + fe ees 
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OB 2 


2 
-$y-4ch 3s eer B25 [riser 
5 “00 
So far everything is exact; Eqs. (2.5), (2.6) and 
(2.7) together with tne first and fourth of 
Eqs. (2.4) are equivalent to the original Maxwell 
Note that in Eqs. (2.6) and (2.7) the 


cancellation of the two large radiation terms has 


foe] 


equations. 


taken place. 
The high-frequency approximation consists now 
in taking the awkward set of equations (2.5), (2.6) 
and (2.7) and ignoring terms in 3/dr, as compared 
with those in d/cdr of the same quantity. 


for example, the first two terms on the right-hand 


side of Eq. (2.5). We see that 
T T 
(2.8) rE = 3 rE cdr >> 2 rE cdr 
9 cdr ) or ) : 
=00 -00 


Consider, 


and so on. We thus obtain 


(rE, ) ~ 

: 1 
wx * ano{ rE,) S SPREE S Sg 
(2.9 


O(rE_) = 
") : _i__ 3E 
5.7 * eno( rE.) = anrdy + sino a 


where we have written 


T 


dE 
(2.10) B def E. - [ a edt 


-00 
for short. We do not drop the smaller integral in 
this case because it introduces no complication. On 
the other hand, it does extend the range in time over 
which the approximation is valid. All transverse 
field components can readily be obtained from 
Eqs. (2.9) and (2.10) once the longitudinal field 
components are known. For these components we have 


the equations 


OE 


pa 1 
c a + 4noE = -knd + 


Yr sin * 


F a 3B, 
(2.113 36 (sin 3 By) ae gt oe 


OB 1 dE, 3 


dos ry = sin 3 E 
c Or ~ Yr sind) oD 56 Teg . 


How these equations are solved depends on whether we 
are concerned with a ground burst or an air burst. 
For the ground burst, the operator, 3/09, becomes 
very large near the ground and this allows an approx- 
imation scheme discussed in detail in Ref. 2. 

In this report we confine our attention to a 
low-altitude air burst, in which case there are 
three possible asymmetries: 

a. Asymmetry of the bomb itself. 


b. Asymmetry resulting from vertical 
gradients of atmospheric properties 
such as density and water vapor 


content. 


c. The geomagnetic field. 
All three of these are small effects, the second be- 
cause atmospheric gradients are small, and the third 
because Compton electron range at low altitude is 
much shorter than its Larmor radius in the geomag- 
netic field. 


are the transverse components. 


As the asymmetries are small, so also 
Moreover the angular 
derivatives are small operators, of order 1/r or 
smaller. The result of tnis is that the transverse 
fields may be neglected in the radial E-equation 


giving 


1 o£. 


(2.12) Sym + nok, = -bnd , 


whicn is immediately integrable by quadratures. 
Foraal integration of the radial B-equation yields 


= 
1 3 
(2.13) rB. = rt oo i E,cdr 
-cO 


Using this in the remaining equations, we see that 
the new terms involving integrals may be dropped for 


exactly the same reason as could the others and we 


have 
d(rE ; 1 3E 
oot 2no( rE,) = -2nrdJ, + 5 39 
(2.14) 
o(rE,) 1 dE 


+ ano{ rE) = -2nrd_ + 


or o 2 sind @° 


Thus our procedure is first to solve Eq. (2.12) for 
En next to calculate E from Eq. (2.10) and then 
Eqs. (2.14) are soluble by quadratures. 

Note that the geomagnetic field enters solely 
through the quantities, Jg and J_, whereas the other 
asymmetries enter solely through 0E/de and dE/d9. 
Furthermore note that E will be a function of atmos- 
phere density, p, of water vapor content, W, and of 
bomb asymmetry, 2. Thus we may write 


dE ap OF , OW OE , ON OF 


(2-15) 39" 36 3p * 30 OW * 30 38 * 


In other words, the right-hand members of Eqs, (2.14) 
can be broken into a series of terms each of which 
expresses a single asymmetry. The same is therefore 
also true of the fields. We shall consider sepa- 
rately the individual asymmetries, knowing that the 
resulting signals may be combined linearly to obtain 
the total signal in the complex real situation. 


III. THE SOURCE FUNCTIONS 


The gamma rays from an explosion result from 
inelastic scattering of neutrons in the bomb materi- 
als themselves and in the air immediately surround- 
ing the explosion. Other sources of gamma radiation 
such as neutron capture and fission fragment decay 
are of too low intensity to be of any importance to 
During the reac- 
tion the neutron population in the bomb rises as eT 
and, after the peak which we conveniently define to 


be at tT = 0, falls exponentially as e5™_ noth a and 


the prompt electromagnetic signal. 


& are of order 10 /sec. The prompt gamma rays follow 


this same time history at the source. Fast neutrons 
which escape scatter in the air and, with every ine- 
It is easily 
seen that this population also decreases exponen- 

tially with time, say as e. 
is much lower than that of bomb materials Ko is much 
less than &, typically around 10°/sec. The air ine- 


-KoT 
lastic gamma rays also vary with time as e . . The 


lastic collision emit a gamma photon. 


T 
» but as air density 


Compton electron current is proportional to the 


gamma-ray flux and therefore varies as 


(3.1) J= -9(r,8,p) f(r) , 


where 


ou for 1r<0O, 
(3.2) tr) = 
-3T “KT 
e + €e for T>0O 


3 


near the explosion. 

The electromagnetic signal, however, comes 
mainly from a region which is 5 to 10 gamma-ray 
mean-free-paths, i.e. a distance of 1 to 2 kilome- 
ters, from the explosion. At such distances scat- 
tering strongly modifies the time history and the 
intensity of the gamma rays and the Compton current. 
A convenient way to describe these modifications is 
in terms of a build up factor, B(r), which is the 
ratio of scattered to direct radiation at a dis- 
tance, r, and of a response function, u(r,t), which 
describes the time history of arrival of scattered 
gamma rays at a distance r from a 5-function source. 


Then at a distance r from the explosion, we have 


(5.3) J = -9(r,*,0)4 f(r) + B(r) 


f(s)u(r- on. 


-0O 


where 9 includes the exponential absorption and 
1pr® attenuation of the unscattered beam. 
A fair approximation to the response function, 


u, is the simple exponential, re » where «, is 


(e) (e) 
constant in +t and varies slowly with r. In this 
case Eq. (2.3) readily works out to be 

«,B 
1+ a ae e for T< 0 
att 
0 
«,B “KAT 
oor 1- g + «Be 
8- Ko 0 
(3.4) J= -9 
1 1 € | 
a+ Ko 8 - Ko Ko 7 Ko 
«, Be -K,T 
+ : om e . fort>0O. 
0” “e 


Typically Ko is of order 10! so that it is consider- 


ably larger than « At the large distances which 


interest us, the mer up factor, B, is quite large, 
of the order 10, and therefore we see that J falis 
as Pee rather then e°7 throughout most of the 
prompt period. 

At the large distances of interest to us, 


electron-ion recombination is completely negligible 


and the conductivity, g, can be separated into 
electronic conductivity, Ce and ionic conductivity, 


G35 governed by the equations 


do, 
ae eh 
(5.5) 4 do, 
i 2 
cd 08g) PAS 


where 6 is the electron attachment rate in air, 6 = 
10°/sec at S.T.P., and A, is a factor depending on 
electron mobility, He» the number, v, of secondary 
electrons made by a primary, and the Compton elec- 
tron patn length in air, f, namely 

(5.6) A vu /t . 

The quantity, Ay» is the same thing but formed with 
twice the ion mobility, 2u,, instead of He (two ion 
species), and a is the ion-ion recombination coef- 
ficient expressed in appropriate units. On our 
time scale it is quite accurate to set a= 0. Be- 
cause the electron attachment time, oc is so short, 
the conductivity, Og» follows the prompt source, a 
-A.d; except for a short period of duration about Bp 
after the peak. Thus o also rises as eT for T< 0 
“SOT at first. Later the conductivity 
levels off to a broad shoulder falling, let us say, 


“KoT 
as e . 


and falls as e 


This leveling off may occur because air 
inelastic gamma rays have caused J to level off. In 
the absence of air inelastic gamma rays, o levels 
off because of the build up of ionic conductivity. 


KoT 


As we indicated above, Koe is only a fair 


approximation to the scattering response function, 
u(r,T). It can be made a good approximation, how- 


ever, if we replace the constant, « by a time 


0’ 
It follows, then, that 


except for a very short period of order, (0<1<1/8), 


dependent function, «(r,rT). 


the conductivity can be accurately described by the 


expressions 
c 
o= 2) (r,9,9)e"" for r< 0 
(3.7) : 
-f nar 
) “Kat 
o = 4s ,(r,8,9) Ke + me for tT>0O 


with « a slowly varying function of r and tT, n a 
slowly varying function of r, and Kp a constant. 
The t > O part of such a curve of o vs. Tt is shown 
in Fig. 1. 
to differ; in fact they must be so chosen that J is 


The quantities, o_ and o,, are allowed 


continuous at + = 0. From Eqs. (3.9) we now find 


(5.8) 


J. =z=- =e -« + «/e)ee ? 


~ “Kot 
+ (Bp - ko) ne ror T>O. 
Continuity at rt = 0 yields 


B-K, 
B+aQ ¥ 


Ko(6 =, * Ko/K) 


p+a oa 


(3.9) ey 


wnere dots mean d/dT and zero subscripts mean the 


value at t= 0. We have defined 


T 


(3.10) wt fase. 


(e) 


when there are many air inelastic gamma rays, the 
electronic conductivity levels off before ionic con- 
ductivity becomes important and we have 

(5.11) = H (with air inelastic) . 


In the absence of air inelastic gamma rays 7 vanishes 


and 
- 


A 
: ie ¥ - 
(5.12) rts f (B - « + &/«)(Ke Mar) 
ie) 


Ay a -p 
= HH [8 - 6) - (B- e*] . 
Thus, as (6 - K)A,/eA, is very small 
ay 
(3.135) qe i, (By - Ky) 1 fy 0 (no air inelastic) 


in the absence of air inelastic gamma rays. The peak 


to shoulder ratio, Ko/n» is normally of the order of 
100 at distances of interest to us. 

With the above general expressions for o and 
for J we can calculate the radial E-field. For this 
we define two auxiliary quantities, S and Es? defined 


by 
es 
(3.14) a kxcodr 
ie) 
and 
def 


(3.15) E, 7 -J/lo . 


From these definitions we readily see that 


(3.16) 8 * nee 2 [a -e"+ 2 (i - oe) | 
2 


fort>O, 


(S is undefined for t < 0), and 


E,= B= (B+ a)/A, for r< 0 


“Kot 

, | (Beet k/eje ¥+(B- Ko)(F/n)e 

(3.17) E_ = a 
- e e + (n/e)e “27 


fortr>O. 


All of our equations so far are accurate except for 
a period of order (1/8) after the peak. During this 
period, that portion of the signal whicn arises from 
asymmetry of E. is very small and sizeable percent- 
age errors here simply do not matter. 
We can now solve Eq. (2.12) for E,- In fact, 

from our definition above of E. we see immediately 
that x! 


= 
eg! uncodr’ 7 a uncodr” 
(3.18) E. =e pa Ee uncodr’ . 


-0O 


As E. = Ey = constant for Tt < O we have immediately 
(3.19) Er) = E,[ - cael fort<0O. 


For t > O we break the range of integration at 1+ = 0 
and we split S into two terms, writing 


(@) 1 
tJ 4 SS 
(5.20) 85 MOL 5 
“KAT 
S, = 4ne 3 {ew- 2. . } 
L « 
2 
and obtain 
. -S(r) 
(5.21) E(r) = E(-o)e 
7 
S)(1) -8)(1’) 
re E,(t‘)e anodr’ , 
fe) 


where E,{-0) is given by evaluating Eq. (3.19) at T = 
Oo. 

Note that Sy becomes much smaller than So after 
a period of a few times (<)7. During this short 
period when the two are comparable, E, varies but 
little and this circumstance enables us to evaluate 
the integral of Eq. (3.21) to good precision, ob- 


taining simply 


(3.22) E(r) = E(-o)eS + g,[2 2 eS ] 


where Eo is the value of E. obtained by evaluating 
Eq. (3.17) at tT = +0. The details of evaluating the 
integral to obtain Eq. (3.22) are all worked out in 
Even Eq. (3.22) can be simplified. As 
we shall see, at times in the neighborhood of the 


Appendix A. 


peak the signal comes from a very narrow region about 
the point r = R, where 2no = 1/a, very nearly. In 
this neighborhood, therefore, #nco/a is very small 
and 


2c 


so that the term in E{-0) may be dropped, giving 
Ss 
es 


As we shall see, carrying the term in E,(0) would 


(3.24) EL = E(1 -e 


introduce no additional complication, but would not 
contribute noticeably to the signal. 

In order to evaluate the integrals required to 
Calculate the radiated fields, it is convenient here 
to define a new function, k, by the relation 


“KAT 


Leh 2 
(3.2) & 98% s/oetne ee ees ; 
neh + ne <2" 
fortT>0O. 


so that Eq. (3.24) may be written as 
: . -ko 
(5.2b) B(r) = E(2 -e)fort>o. 


We see that k depends only weakly on position but 
strongly on T. 

To complete our evaluation of tne source func- 
tion we need the time integral of Eq. (5.26). From 
kq.(3.16) we see that S, i.e. ko, is essentially 
constant in +t except for the first cuuple of yener- 


ations after the peak. Therefore we can write 


T 


(3.27) “Ky 6 


E.cdt = E(l-e T 
r (e) 


to good precision, for when tT is so small that ko is 
varying this expression does not differ significant- 
ly from zero. Thus we have 
T 
dE 


def r 
or edt 


(3.28) E = ae 


fe) 
OE 
- 7 ¢) -ko -ko 
= (= cT  }o -e )+ E_Zetkoe $ 


where Z is defined by 


zy def 1 o{ke) 
ko or ~ 


Clearly Z is about 1/\ and the correction term to z. 


(3.29) 


would not be needed for ct << x. We keep the correc- 
tion term in order not to be restricted by such a 

severe limitation. Tne source function for atmos- 
pheric and for bomb asymmetries is obtained by dif- 


ferentiating Eq. (3.28), and is 


2 
= oE o°E 
OE _ fe) 0) ko 
(9220) 39 E et sas |[2 = 
OE dE 
(e] 0 fo¥A -ko 
+ XE - er(n x - 2 30 Eo 22) koe 


+ Ae zen] ke - (xo)*]e*9 : 


where we have written for snort 


(5.51) 9 S8f 2 Sk) 

The terms proportional to T, which arise from 

Eq. (3.27), we shall call the secular terms; they 
are clearly of importance only when c1/a is appre- 
ciable. 

In deriving Eq. (3.30) we have allowed specifi- 
cally for all effects of gamma-ray scattering except 
for one, namely the change of angular distribution 
of Compton electrons with distarce and with time. 
This effect can, however, be included by simply re- 
defining the quantity, v, the effective number of 
secondary electrons per Compton electron, to be a 
function of t and of r. This introduces additional 
rt and r dependence of Es) but we already have specif- 
ically allowed E. to depend on these quantities. 
Thus, with proper interpretation, our results are 


quite general. 
IV. SPECIFIC SOLUTIONS 


We have seen that the geomagnetic signal can be 
discussed separately from the other signals. As 
this signal has been discussed quite well elsewhere 


(5,4,5) we content ourselves 


and is well understood, 
here with deriving the basic formula which we do 
simply for completeness. We consider a perfectly 
symmetrical explosion in a perfectly uniform atmos- 
phere but in the presence of a uniform magnetic 
field, Bo: 
or ~By we find 
we 


Choosing our polar axis in the direction 


se 


(4.1) Jo 3g J sins, J,=0, 


) 


wnere £ is Compton electron range, a is its giro 
radius in the geomagnetic field and J is the radial 


Compton current. Our only surviving field equation 


is 
o(rE,) ! 
(4.2) ser 2no( rk.) = -anrd (32) sins, 
whose solution is 
r 
(4.5) rE > 4 sin’ af rE,e “2nodr ; 


fe) 


At large distances we can set r = oo in the rignt- 


hand side and obtain for the radiated signal 


co 
£ sin’ - 
(4.5) TED = =a ‘i rE.e Xonodr . 


fe) 


The factor, anze*, has a very sharp maximum at r = 
Ry the point at which 

1 do 
(4.6) ergs Ze -sy- 
This fact enables us to evaluate tne integral. The 
details are given in Appendix B and, to lowest order 
in AR,, we obtain tne well-known formula 


£ sin 9 
(4.7) TED = — ta 


* RE,(R.,7) ‘ 

So long as R, is considerably larger than A, Eq.(4.6) 
above indicates that it is proportional to log o. 

E, is constant during the Q-phase and changes rapid- 
ly at the peak, dropping to a lower nearly constant 
value during the x-phase. Thus the geomagnetic sig- 
nal is a short, sharp spike having its peak slightly 
before o reacnes maximum. 

We now address ourselves in somewhat more detail 
to the other asymmetries. We may ignore geomagnetic 
effects, setting Jy = Jy = 0. It is convenient to 
choose our polar axis vertical. If then we suppose 
tnat the bomb itself is azimuthally symmetric, so 


will be the field equations, which reduce to 


(rE, ) 
(4.8) ase + 2no(rE,) = zo ; 


which has the solution 


= e& 1 3E .-x 
(4.9) rE, = e [3d dr , 


X being the quantity defined above in Eq. (4.4). We 
confine our attention to the distant field, writing 


and shall refer to this as the radiated signal. If 
we now substitute for dE/88 the expression given by 
Eq. (5.30) we obtain 


co 
(4.11) rE, = f @ [> - oe] ear 
1 
0 
Cc co 
“f opkoe KO *ar + f | ko - (x0)? Je"KOKar- , 
0 0 
where Py> P55 Ds are given by 
[Eo a°E 
°1 > B/3e ~ °T dre |’ 
dE dE 
ah 3 0 ) 
Pe Gees a, + er( 2 SE + tae - 9 3e) , 
1 
5 = 3 NE, zZer . 


We have defined Z and 9 in Eqs, (3.29) and (3.31). 

Now consider for a moment the second term of 
Eq. (4.11). The factor, ka exp [-ks - x],is a rapidly 
varying function of r with a sharp maximum at r = Ry 


that value of r for which 


def 


(4.l2a) Oncy = Z, y = 1+kZ/en. 


On the other hand Po is a slowly varying function of 
r and this enables one to evaluate the integral with- 
out specifying the functional forms of » and of o, 
In Appendix B this is worked out in detail, as are 
also the first and third integrals. 


expansions are given out to third order in (a/R) - 


In this Appendix 


As this quantity is normally about 1/lo, it normally 
suffices to take tne first terms of these expansions 


obtaining the result 


9k Qk 
(4.13) rt, - [pio y + + | 7 
reR 


If higher order accuracy is required, additional cor- 
rection terms will be found in Appendix B, along with 


error estimates. 


The quantity. k, [Eq. (3.20)] is a rapidly vary- 
ing function of t and, at most, a Slowly varying 
function of r. 
most linearly with t. Normally ?, is of the samc 


The m’s,on the other hand, vary at 


order as the secular term of Po and, as y becomes 
quite large, the term in ?, is a relatively small 
correction. 

We can now dispose of the Q-phase in 4 general 
way. During this period we have k = 4nc/a@ and Z is 


very nearly 1/A so that 


2c ec/a 
TE, = NP, log (: +H ) +e, 1 + Sc/on 


+o ec/a 


? (1 + 2c /an)” , 


(4.18) 


which 
) 


are proportional to Tt are completely negligible at 


But 2c/oy << 1 and the terms of 1» Pp. ® 
this stage. Thus we may expand Eq. (4.14) and 


obtain 


(4.15) ig Sy{k tees fort <0, 
where E, is the (constant) value of E, when + < 0. 
The signals which come later are Clearly of order 
Ep which is much larger than 2cE,/a. We can there- 
fore completely ignore the Q-stage signal. 

For t > 0 the radiated signal behaves essen- 
tially like the first two terms of Eq. (4.13), the 
third remaining always a small correction. Thus, 
aside from the secular terms, which are important 
only when Tt approaches a microsecond or so, the 
signal is given by a linear combination of the two 
functions of +, log y and Zk/eny. Over most of the 
range of T, zk/2x is large and these become slowly 
varying functions of k and hence insensitive to its 
details. 


quantities not as functions of real time, Tt, but 


This becomes clearest if we examine these 


rather as functions of generation number, [, which 


we define as 
aff sien of ] 
(4.16) © ' sign 7, 108 (o,cax/0(7) rR,” 


In Fig. 2 we have plotted the two functions 


F(P) = log (1 + k/2na) 


k/2nr 
G(r) = 14 jexx 


as functions of the generation number, I. 


(4.17) 


Each was 


calculated for two different cases, namely: 


a. The solid curves are computed from 
the o vs. Tt history of Mg. 1. This 
case is described by a prompt burst 
with « varying from an initial value 
of 2.8 x io! to a final value of 
5.8 x 10°, plus a long tail of air 
inelastic gammas for which Ko = 
ux lo’. The peak to shoulder ratio, 


(xo/n), is 70. 


b. The dotted curves are computed for a 
prompt burst of constant kx, k # k& > 
2.24 x 10" and without any long tail 


to o. 


The curves of case a. correspond to a o history 
which is far from being a simple exponential, yet 
the dotted and the solid curves lie remarkably close 
to one another. It is clear that for many purposes 
“universal” curves for the two functions, F and G, 
would suffice. In addition, of course, one would 
need ag vs. tT curve to translate generation number, 
[., into real time and for times not much less than 
A/e the secular terms in P) and Po would have to be 
included. 

We shall now proceed to discuss in detail 
special cases. We have in every case retained all 
the secular terms, including the small one in , but 
have calculated only to lowest order in (A/R,), 
which we have supposed to be 1/10. In all cases 
calculations were made for the two o vs. rt histories 
designated as a. and b. above. Curves are shown 
(case a-solid, case b-dotted) for normalized signal 
strength vs. generation number [ and also vs. real 
time, tT. For the quantity, Z, itself, we have used 
simply 1/A, but have used higher order terms to com- 
pute 3Z/de. 


A. Bomb Asymmetry Alone 
lf the only asymmetry is that of the bomb it- 


self, then J and o have identical angular dependences, 


LO 


which may be written as 
J= r,t) > of@) 


o = o(r,t) - ae) 
Q = dw/ade . 


(4.18) 


Thus Ey is independent of 8; as there is no atmos- 
pheric gradient it is also independent of r and 
Eq. (4.13) reduces to 


1 k/2mr 
(4.19) rEg = 5 EQ > ax {1+ + fix} - 


In case b. the secular term is completely unimpor- 
tant, in case a. it never exceeds 0.1, and this 
only after 1 microsecond. Thus this signal is 
essentially described by the function, G. Curves 


are shown in Figs. 3 and 4 giving rE, vs. gener- 


8 
ation number, [, and real time, Tr, respectively. 


B. Gradient in Water Vapor Concentration 

Next consider a perfectly symmetric explosion 
in an atmosphere of uniform density. However, owing 
to a gradient in water vapor concentration, the 
electron mobility, and hence o, varies with altitude. 
If, for convenience, we take the variation to be of 


the form 

Pa 
(4.20) He = Ho exXP {-§ cos 8} 
we find 


R, sin 8 an R. cos 4 
O08) © 2 35 Fo E(t A )| 


R_ sin 8 
Ss 2a, 2 cos 8 
~,(R,) -— a Eo - (+28 i ees 


R. sin 9 ae 
95 (BR) 83a Bg [isd 


(4.21) 


Assuming one is near enough to the equatorial plane 
that R. cos S/n << 1, the radiated signal is 


R. sin 9 
(4.22) rE, = AE)| sp 


{p - Aen See -(.- =) log (1+ Ken) 


to lowest order in APR. Here the secular term is 


of order ct/X and quite important at 1 microsecond. 


Curves of TE, 


time, T, are shown in Figs. 5 and 6, respectively. 


vs. generation number, [, and real 


C. Atmospheric Density Gradient 
At low altitudes air density varies exponen- 
tially with the altitude so that we may write 


(4.25) p= Po eXP {-r cos a/u } 5 


The quantity, H, is called the relaxation height and 
may be taken to be about 9 kilometers. In this sub- 
section we shall always use the subscript O to indi- 
cate the value of a quantity at burst altitude. As 
electron attachment is a three-body process, the 


attachment coefficient, B, varies asp so that 


(4.24) -2r cos $/H . 


B = Bye 
Now writing 


(4.25) J = -9(r,8) - Fir) , 


we can write for 9, 


b 
(20) 9+ BE 58] ow [-E(2- 4) 
0 ! 


So long as we confine ourselves to regions suffi- 
ciently near the equatorial plane that R, cos e/H << 
1. The bomb yield expressed in appropriate units is 
Y and b is a constant arising from an empirical fit 
to the build-up factor. 
Appendix C. 


The details are given in 


In calculating derivatives of BE. from Eq. (3.15) 
we shall ignore any dependence of x, n, and A, on 


Position as we do not at present know these. Thus 


we write 
(4.27) OF, _ 1 a6, Po. ar sins 
: 3e A 8 H 7 
e e 
oF. = ._0 2 cos @ 
Or H fe 


with Eqs. (4.26) and (4.27) and the relation, o = 
J/E,; we can calculate the functions P > Pp, and 3: 
When the result of this calculation is set into 

Eq. (4.13) we obtain 


er/n, 


R k/2nr, 
( i+ k/2nr, 


s Ay 4Bo 
XT ¥ kee, aida ey Wace Oa 


8 
etfo, 4B k/2nr) ik fs 5a 
2 So 7, 2 Aa tT ? 
R, AE, | 1+ Kenn, Ro| ALE. 


correct through first order in A/R,. As Bo/ALE, al 
and \/R. » 1/10 we see that the secular terms are 
again of order ct/A. Curves of this signal vs. 
generation number, [, and real time, T, are shown in 


Figs. 7 and 8, respectively. 
V. DISCUSSION AND CONCLUSIONS 


We have shown how the early portion of the radio 
flash from a low-altitude air burst may be computed. 
Use of the high-frequency approximation enables us 
to obtain simple analytic expressions for the signal 
starting from a quite general o vs. Tt history. A 
question still unanswered remains, namely, for how 
long a time is the approximation valid? To answer 
this question we have used the solutions found to 
estimate those terms of Maxwell's equations which 
were originally dropped. Inserting these estimates 
into Maxwell’s equations, one can solve and obtain 
first-order corrections to the high-frequency ap- 
proximation. This process is carried out in detail 
in Appendix D, where it is shown that the fractional 
error involved in using the high-frequency approxi- 
mations is about cr/R., R. being our usual signal 
radius defined by Eq. (4.12a), tnat is to say, the 
approximation is good so long as 
(5.1) er<< R.- 

Typically this means that we can have confidence in 
about the first microsecond of the computed wave- 
form. It should be pointed out here that condition 
(5.1) was achieved by including the secular terms, 
which our earlier theory did not do. Without these 


secular terms the much more stringent condition 


12 


cr << \ must be fulfilled. 

The second conclusion we wish to point out is 
the limited variety of possible pulse shapes. 
Leaving aside the geomagnetic signal, we saw that 
all other asymmetries lead roughly to signals which 
are linear combinations of two basic pulse shapes, 
described by the functions, F and G, of kq. (4.17). 
Because of the presence of the secular terms, this 
is not strictly true; nevertheless, it is true ina 
practical sense as illustrated in Fig. 9. In this 
figure we have plotted as a solid line the atmos- 
pheric asymmetry signal of Fig. 8 (solid line), 
together with the dotted curve which is a linear 
combination of the bomb asymmetry signal and water 
vapor gradient signal (solid lines) of Figs. 4 and 
9%. The linear combination chosen is an improbable 
one, but it does show that there are essentially 
only two independent pulse shapes. This should, of 
course, have been obvious from the beginning, for 
there are only two basic asymmetries, that of J and 
that of o. Practically, this means that inclusion 
of new effects, such as the spatial dependence of 
A, and of «, will introduce no new pulse shapes. It 


also means that it is impossible to sort out the var- 


ious asymmetries from a given signal. 
There is one effect which we have not included 


in our calculations which might lead to an additional 


wiggle in the curve, and tnat is tne effect of dif- 
ferential scattering. Owing to the density gradient 
of the atmosphere more of the gamma rays arriving at 
some distant point r have been scattered downward 
than upward. This will produce a net current of 
Compton electrons downward at the point in question, 
i.e.,a negative Jo: The signal produced by this Jo 
is calculated exactly as was the geomagnetic signal 
and will in fact resemble the latter signal, with 
two important differences: 


a. Because only scattered gammas contribute, 
the signal will be delayed and broadened 
as compared with the geomagnetic signal. 


b. The polarity of the signal will be posi- 
tive regardless of magnetic bearing. 


A rough estimate indicates that the amplitude of 
this signal is small, but further consideration of 
it is surely viepanvea le” 


Noniintar effects have been omitted frum our 
theory; in fact they are neyligibly small. As the 
siynal comcy from a thin shell centered on r = Ko 
it suffices tu confine our attention to this regiun. 
Using Chr furmulas we have developed for EL. at tnis 
position we see that this quantity is simply not 
larye mnough to extract a significant fraction of 
the Compton clectrun energy. The transverse fields 
are ineffective in producing a transverse J for tne 
simple reason that Ey» a are very nearly a pure 
radiation field so that their effects on the Compton 
electrons very nearly cancel. The only nonlinear 
effect remaining is the field dependence of the 
electron mobility, He» and the attachment coefficient, 
B. Over a wide range of field strength, B and He 
depend on the electric field as ( 


(5.2) bos HWE, Bx BAWE, 


where numerically Ho = 10° and By = 108 when E is 
expressed in e.s.u. During the x-phase, the elec- 
tric conductivity depends on these parameters as 


; is 
(9.5) gee = ————.. 


It turns out that « is so much smaller than B that 
the nonlinearity of He is essentially cancelled by 
that of B and the field dependence of o is s0 weak 
that it can be safely ignored. During the a@-phase 
we do not necessarily have a@ << B. On the other 
hand the radial E-field at r = R. is very small 
during this phase, around 0.03 e.s.u. or less. At 
such low field strength the conduction electrons are 
essentially thermal and He» B are no longer field 
dependent. 


APPENDIX A - THE RADIAL E-FIELD 
According to Eq. (3.21) the problem of calcu- 


lating the radial E-field is solved once we evalu- 
ate the integral 


co fee] 


(8.15) ze f pe Vanodr = f (3) ewer) 7 


r r 
which partially integrates to yield 


co 


(B.16) ita ehe f (8) (rise yar: 


r 


Using Eq. (B.11), this may be written as 


(B.17) 1s e isso tEs 


—. 
8 
b 
nl+ 
< 
on~-~ 
«138 
~~". 


“Vv 
l-e ’ 
Patel 


and we can integrate by parts again obtaining 


(B.18) I -2(1 a e¥) + =P) ty + log y + E,(v)] 


: {Heo Ly + logy + E(w) }ar , 


where y is the Euler-Mascheroni constant, 7 = 
0.5772156 °°°. 
Now put r = R, into Eqs. (B.13) and (B.18) and 


add. We obtain 
cO 
-¥ Q Lee y’ 
B.1 e 'Onodr ={—<1 + —~>-=— 
(B.19) f° [24 ZAG ) 
ie) 


(7 + 108 v}| Hee 
rsk 


8 
The remainder, indicated above by the dots, is readily 
estimated with the aid of the integrals of Kqx. (B.t5) 


and (B.18). The outer region contributes 


¥ 
\/ S00 n-l_ n-l 
1+ w/fl + wi[? So y dy 
oad Zz z ly nnf 
rhe n=1 
id) 


co 
1+ ¢ + Ww el’) . S (-)2 V8 
Zz z y 2 
r=R * is 
s 


Similarly the inner region contributes 


s é 
E,(¥) 
a 1+ w/l + wo)’ 1 ‘ 
(B21) €sn4 = f z (4 2/' woe 
6 
foe) 
‘ 
Lew Lee] B,(v) di 
31 Zz z Ly v ; 
r=R 


s 
s 


Chandrasekhar’ ©? gives as tne value of this integral 


co 
er ee Ie a 
i. n=1 


(B.22) 


whence we find 


2 ‘ 
2 is Py gle Cee) } 
(B23) Emoe Ez IB i 


very nearly. 


Zz 


The error is worst at late times such 
that kz/en is large. 


error is about 


Setting z y 1/A the fractional 


pio «2 {- HE), Ip 


s 


(B25) fag 


tur functions, op, which vary as some low power ol r 
this error is about lof if K, is as small as 5). 


a) 


For yields of reasonable magnitude R, is more like 

10X and the error becomes quite small. If to 

Eq. (B.19) we add the correction term, Eq. (B.25), 

the remaining error is of order (AR)? and thus 

small even for such tiny yields as R, = 3dr. 
Applying Eq. (B.35) to Eq. (B.19) yields 


similarly, from Eq. (B.4) we have 


a se 
en | . , 
y r=R 
8 


In order to calculate second-order corrections 


8 


( B.25) oen(o)@e Var 


iy tne field quantities. we also need I, and I, as 


2 
functions of a running upper limit; I, however is 
needed only in correction terms, which is why 

it is given above to first order only. We have 
already calculated I(r) and we can see that, as 
g(r) is presumed to vary slowly with r, we have to 


first order 


or) _-¥(r) 
Be e for rs R. 


r 


(3.26) [sear 


? 9 


ac 
Ys 


er) 


for r2R_, 
s 


where subscript s means evaluation at r = Ro: Now 
frum Eq. (8.5) it follows that 


 [a(22) - a 


(8.27) J Uli ee ae 


(e) 


Lu 


As these quantities are used only in correction 


terms,, we do not need higher order accuracy. 


APPENDIX C - CURRENT DENSITY IN A NONUNIFORM ATMOS- 
PHERE 


If we write out Eq. (3.3) specifically for a 


uniform atmosphere we should have 


fc 1) Je- alin {£(1) + Br)F(r)} 
. 2 B 
r 


where we have written F(t) for the faltung integral 
of Eq. (3.3) and Y represents the yield in appro- 
priate units; A is of course the gamma ray mean- 
free-path. ‘The build-up factor is quite weld. repre- 
sented empirically by 


r Me. 
(c.2) Kr) = oF) , 
with b a constant. For distances such as r « 10d, 
it turns ont that B >> 1 and we can in fact write: 


(C.3) J= -9(r) > Nr), 


’ ves 


where 


(ony ge Mery 4 ok 


Sometimes in fitting Kr), an exponential factor is 
included on the right-hand side of '£q. {C.2), but 
this merely leads to a redefinition of A in our 
Eq. (C.4). ; 

In an exponential atmosphere, » is no longer a 
constant but is given by 


os 


a) 


er cos 3/H Fe * 


(C.2) $= 


In Eq. (C.4) r/x means the distance measured in mean- 
free-paths and tnerefore, in case of a variable den- 


sity we must make the replacement 


a r dr 
(C.6) = — 


- H -r cos 8/H 
x Ay cos 8 {l-e J 


(e) 


When this replacement is made, Eq. (C.4) is found to 
correspond with 


You 


r A, cos $ 


ie) 
H -r cos 6/H 
op [- Bogie | 


in the case of the atmosphere of variable density. 


(C.7) 9= (1-e% 78 8/4 


In the end this quantity and its derivatives 
will be evaluated at r = R.- At low altitudes R, is 
considerably less than H. Moreover the signal is 
frequently observed near the equatorial plane so 


that 
(c.8) R, cos @ << H . 


When this is true we can expand some of the exponen- 
tials simplifying 9 to 


(c.9) 9 = ral - 298 * Jew] - Z(1- 


which is Eq. (4.26) of the main text. 


APPENDIX D - LIMIT OF VALIDITY OF THE HIGH-FREQUENCY 
APPROXIMATION 


In making the hign-frequency approximation, cer- 
tain terms were dropped from Maxwell's equations be- 
Now that 
we have obtained such a zero-order approximation to 


cause they are unimportant at early times. 


the fields, we can use these results to calculate tne 
values of the terms dropped from Maxwell’s equations. 
These additional source terms enable us to calculate 
first-order corrections to the fields, and in princi- 


ple one could continue, writing 


= (0) + (1) toe 


~~ 
o 
° 
w 
~~ 
ty 
i) 


bg 
cy 


(0) (1) 

B + 5 a a 
9 Bp 

We shall carry out this program as far as the cor- 


rection gf) to the radiated field. Our aim in 
this calculation is not to achieve higner order 
accuracy of the fields but rather to define the 
time interval for which the zero-order solution is 
valid. Thus we shall content ourselves with some- 
what rough evaluations of the integrals wuich will 
arise. 

Owing to the term proportional to ct, tne term 
in Pp is normally larger than the rest. AS we are 
only estimating an error it suffices to consider 
this term alone, i.e. we may choose a model in 


which 


vy Ps =O, 
(D.2) 
1 
P= 5 ME (1 + et/r) , 


typical of our worst cases. To achieve tnis we nave 
0E,/04 = O, whence from Eqs. (3.20) and (2.9) we ob- 


tain 


1 
(D.3) rBy =} xspa(s + SE) Pa 


IN 
and 
(D.4) = 3 { sin ) | 
r sin 6 
é 
, Ee (2 4 St ) k/enr : ect ke eke 
ore ~/1 + k/enr x ? 
where 
, , def 1 3 Q si 
(D.5) af = 3in 6 30 (Q sine) . 
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In all our cases 9 is of the form € sin @ whence 
i’ = 2€ cos 6, where € is a small asymmetry factor 
whose exact form depends on the type of asymmetry 
under consideration. In all cases, however, 2! << 
1. Equation (D.4) gives the source function for 


(1) 


the correction EL ‘to E rr 2-e- 


ut fit by (sin jek NS cay. F) 


sin $ 


- : 
=~ —2 -ko St) k/2nr- 
ao f{o +S) Fen 
0 


“(pe) rE 


(11) 2 
@ 


plus terms of higher order in (a/R) - 


Partial inte- 


gration and the results of Appendix B show that the 


above integral is simply 1/R, to lowest order and we 


have 
rE, 22 
wr (11) | “o 
(D.12) rE, R, cr += |, 
whereas 
eer koe edt 
ae . AE,2 
: “(0) ) = | k/ena 
(D.13) TES ~-3[{1+¢ ea - 


Now the quantity k/@nay is Always Yess than unity 


but is very nearly one for all appreciable values 


€learly the correction term is 


negligible for small 


values of T; it remains negligible just so long as 


of +t. Similarly, ko is nearly constant Tor’ such’ 
values of + and we have ee 
tc 
; (D.14) 3a R, << 1 
XE nf 
(D.7) e}) — onl cr + = (1 + axe) | ; 
er ; (11) 


to adequate accuracy, and tnerefore 


In addition to ES 


order correction ts a2 


there is another first- 
, which has as its source 


the term 
(1) 
3k E,W 
(D.8) 5 36 = = etl er + (2 + 2a] (2) 5 T 
= (D.15) > = 33, rE() car : 
or 
+ o(e*) 


The terms in o( €°) vanish identically when 9 = 


and are in general small. 


n/2 
We have defined 
(D.9) sh = dst/de . 


This correction to E. induces a correction to 
Eg which we shall call EA 1). 


; (1) 
dE 
(D.10) re(1}) -f = e “ar . 


betting Eq. (D.8) into (D.10) gives 


It is given by 


This is readily calculated to be 


ey rE, 22 
(D.1b) >. “2 (ce + = 


2 
a2 3( za 
{25+ an 200 


(25 4 
© 


“ef2) = S 


(D.17) 


Jers 
2a ko - oats} 
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When we set Eq. (D.16) into Eq. (D.17), all of the 
integrals are given in Appendix B except for that 


involving the term in npr. 


we have 


Writing y = ko + X, 


(D.18) 


as the result of partial integration. 
Eq. (B.19) of Appendix B, we obtain 


Now using 


1 


(D.19) eYar = yre , 


when we make use or the fact that V, = 1. Using 
this result plus others from Appendix B we obtain 


2 
XE,2 22 2 
(D.20) re 12) = m- (< ¢ Se ) ye? 2 og 


+ oben dz, 
1 + k/2nr} or 


2a _ 2 kann 


i (1 + k/2na)® 


DIN 


r=R 
s 


The worst case is that of the atmospheric density 
gradient for which Na r, Then the term involving 


zon, 


(v.21) 2982S 
Ss 


dominates the rest and we have 


E 
CT 


() 
(D.22) sa. 
BO) eR. 


From Eqs. (D.14) and (D.22) we can conclude that the 
high-frequency approximation is valid for all times 
such that 


(D.23) ct << R, - 


Note that this conclusion is essentially independent 
of the magnitude of the secular term. 
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Fig. 9. linear dependence of the various signals. 
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